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We investigate the x-ray photoelectron spectroscopy (XPS) of nominally and n-doped dP tran¬ 
sition metal oxides including Nb02, SrVOa, and LaTiOs (nominally d^), as well as n-doped SrTiOa 
(nominally d°). In the case of single phase d^ oxides, we find that the XPS spectra (specifically 
photoelectrons from Nb 3d, V 2p, Ti 2p core levels) all display at least two, and sometimes three 
distinct components, which can be consistently identified as d°, d^, and d^ oxidation states (with 
decreasing order in binding energy). Electron doping increases the d^ component but decreases the 
d^ component, whereas hole doping reverses this trend; a single d^ peak is never observed, and the 
d° peak is always present even in phase-pure samples. In the case of n-doped SrTiOa, the d^ com¬ 
ponent appears as a weak shoulder with respect to the main d° peak. We argue that these multiple 
peaks should be understood as being due to the final-state effect and are intrinsic to the materials. 

Their presence does not necessarily imply the existence of spatially localized ions of different oxi¬ 
dation states nor of separate phases. A simple model is provided to illustrate this interpretation, 
and several experiments are discussed accordingly. The key parameter to determine the relative 
importance between the initial-state and final-state effects is also pointed out. 

PACS numbers: 31.15.A-,71.55.-i,73.20.hb 


I. INTRODUCTION 

X-ray photoelectron spectroscopy (XPS) is a very common in-situ and ex-situ tool used in modern laboratories 
to probe the stoichiometry of a given material, as well as the oxidation states and local chemical environment of a 
given element [ll-Q . As the core levels of different chemical elements easily differ by tens to hundreds of electron-volts 
(eV), the peaks in the photoelectron distribution as a function of kinetic energy provide us with information on 
which chemical elements are present and, to a very good approximation, their relative abundance in a sample. When 
focusing on the photoelectron signals coming from one particular core level of one particular element, the different 
local environments around the targeted ions can result in a multi-peak structure, typically within an energy range 
of about 10 eV, from which the oxidation states of the probed element can be inferred [5|-l7|. A more sophisticated 
aspect of XPS is the electron screening due to the created core hole [^: once a photoelectron is generated, the sample 
is left with a core hole (positively charged) that modifies the potential of valence electrons. The response of valence 
electrons to the core hole is usually referred to as the final-state effect, in the sense that the observed spectrum does 
not really correspond to that of the neutral sample before being irradiated, but rather to the energy spectrum in the 
presence of a core hole. The typical lifetime of a core hole is about 10“^^ s [^, which results in an energy broadening 
of ^0.1 eV. Accordingly, peak features that are larger than 0.1 eV in the core-hole spectrum can, in principle, be 
observed and resolved. 

The final-state effect introduces even more features and complexities to the XPS spectrum, as electron correlation 
is essential to the process of core-hole screening. For example, the XPS spectra of a metallic system typically has 
an asymmetric shape (orthogonality catastrophe) when taking the scattering of the core-hole potential into account 
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-[lOj . In addition, if the targeted ion has degenerate localized orbitals (such as 3d or 4/ orbitals) in a metallic phase, 
a uniform system also displays multiple XPS peaks. To properly describe such systems theoretically, an Anderson 
impurity model including both localized correlated orbitals and uncorrelated bath orbitals is required ^,11, 1^. For 


the transition metal (TM) oxides, the valence states have to include both oxygen 2p and TM d orbitals, as their 
energy difference and their mutual hopping amplitude are comparable in energy. Therefore, a minimal model for XPS 
spectra of transition metal oxides includes a TM-Oe cluster [,‘lL ll,‘ll4l5l |. Although complicated, once the XPS spectrum 
is properly interpreted, it provides a quite good estimate of material-specific parameters such as inter-site hopping 
amplitude t and Hubbard on-site repulsion U. 

In this paper, we reexamine the origin of the multi-peak structure in the XPS spectra of nominally d^ transition 
metal oxides including Nb02, SrVOa 16|, and LaTiOa, as well as that of lightly n-doped SrTiOa (STO) 


17- 


In particular, we propose a cluster-bath model and argue that it is the final-state effect rather than the presence of 
multiple oxidation states that accounts for the observed multi-peak XPS structure in these materials. Based on our 
interpretation, the multiple XPS peaks are intrinsic to the materials, and do not necessarily imply the existence of 
spatially localized ions with different oxidation states or of separate phases. The rest of the paper is organized as 
follows. In Section II we give a brief overview of the XPS core level spectra of these four oxides. In particular, we 
distinguish between the initial-state effect and final-state effect. In Section III we present our experimental results 
and point out their common features and their implications. In Section IV we provide a simple model to illustrate the 
final-state effect, which is crucial to reconciling the seemingly conflicting observations. Several experimental results 
are discussed accordingly. The key dimensionless parameter to determine the relative importance between initial-state 
and final-state effects is identified. A brief conclusion is given in Section V. In the Appendices we provide the details 
of our calculations. 


II. OVERVIEW OF XPS 


In an XPS experiment, photons of energy hv are directed to the sample and photoelectrons of kinetic energy Ekin 
come out [see Fig. [T][a)]. Energy conservation requires that 


hv -\- Eas{N) = Ekin + Ecore {N-l) + cj,. 


( 1 ) 


Here Egs{N) is the ground state energy of the sample with the filled core level, Ecore[N — 1) is the energy with a 
core hole {N — 1 is used to denote the presence of a core hole), and </> is the work function. By shifting the kinetic 
energy by Ekin —>■ u; = Ekin + </> — hv, the photoelectron intensity as a function of w is given by 


p(cc) = ^ \{n{N - l)|c|G'5)p x 6{uj - [Egs{N) - E.^reA^ - 1)]) 

n 

= -(G5|ct [w - {Egs{N) - HtA - idA c\GS). 


( 2 ) 


Here is the creation operator of a core electron, |OS') and \n{N — 1)) are, respectively, the ground state without 
a core hole, and eigenstates with a core hole [Sj, y, ll2|. Once Eltot is specified, the second line of Eq. ([2]) is used 
to compute the XPS spectra. Note that p(w) is non-zero only when w = EgsA) — Ecore{N — 1). What the XPS 
spectrum reflects is the core-hole energy spectrum weighted by the matrix element \{n{N — l)|c|GS')p. The XPS 
spectrum is also routinely plotted as a function of binding energy Eb, defined as Eb = hv — (p — Ekin = — w H- 
Por the purpose of this work, the constant energy shift is not important and we focus only on the dependence of the 
spectrum on the “relative binding energy” or “relative kinetic energy”. 

Conventionally, one distinguishes between the initial-state and final-state effects in the XPS spectrum [^, |^. For 
the initial-state effect [Fig.[I][b)], the valence electrons are not affected by the created core hole. In this case the XPS 
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peak position is determined by the core-level energy tc only. Within this scenario, any observed multi-peak structure 
in the measured XPS spectrum implies that targeted ions (where the photoelectrons are ejected from) experience 
different environments within the same sample. For example at the Si/Si02 interface, the observed multiple peaks in 
the Si 2p spectrum, which corresponds to different Si oxidation states (from Si®’*' to Si^’*'), are used to deduce and 
quantify the formation of SiO^, at the interface Q. For the final-state effect [Fig. [Ijc)], the valence electrons do feel 
and respond to the potential caused by the creation of a core hole. In this case a spatially uniform system can also 
lead to additional peak structure around Cc in the XPS spectrum. A classic example is CeNiy, which is a nominally 
/° material but displays three XPS peaks (from Ce Sd core level), identified as /°, /^, p [2l|. It was realized by 
Kotani and Toyozawa [ll] . , and by Gunnarsson and Schonhammer [l^ that the multiple peaks in this material 

originate from the final-state effect, where the valence electrons response to the presence of a core hole, especially 
the core-hole-induced energy change of Ce 4/ levels, plays an important role. Simply put, for the initial-state effect, 
the ions of different nominal charges preexist in the sample; for the final-state effect, the ions of different nominal 
charges are created after the applying photons produce core holes. We believe the experimentally observed multi-peak 
structure in nominally and n-doped d° transition metal oxides should be understood as being due to the final-state 
effect. In the following we shall provide our experimental and theoretical analysis that leads to this conclusion. The 
key parameter determining the relative importance between initial-state and final state effects will be discussed in 
Section IV.E. 


III. EXPERIMENTS AND KEY FEATURES 


In order to properly analyze the intrinsic XPS spectra of nominally d^ transition metal oxides, we need to be able 
to grow single phase, crystalline layers of these materials and then measure their XPS spectra without exposing the 
samples to air, as these materials are not thermodynamically stable in the ambient and will slowly oxidize. The samples 
of Nb02, SrVOa, and LaTiOa, as well as SrTiOa with several dopants, are grown in a molecular beam epitaxy (MBE) 
chamber and then transferred in situ to a high resolution photoemission chamber. The two chambers are connected 
by an ultrahigh vacuum transfer line with a base pressure of < 1 x 10“® Torr, allowing for sample transfer between 
the growth and analysis chamber within 5 min. The photoemission chamber consists of a monochromated A1 Ka 
photon source {hv = 1486.6 eV) and a VG Scienta R3000 analyzer. XPS spectra of the valence band, O Is, Nb 3d, 
V 2p, Ti 2p, Sr 3d, and La 3d are taken (as appropriate) at a pass energy of 100 eV with an analyzer slit setting of 
0.4 mm, resulting in an overall instrumental resolution of 350 meV (primarily limited by the energy resolution of the 
x-ray source). The analyzer is calibrated such that the Fermi level of a clean silver foil is at a binding energy of 0.00 
eV and the Ag 3^5/2 core level is at 368.28 eV. 

Undoped SrTiOa is nominally d^, while the remaining three materials are nominally d^: in the ionic limit, SrTiOa 
has no electron occupying the Ti 3d orbital; Nb02 has one electron occupying the Nb 4d orbital; SrVOa and LaTiOa 
have one electron occupying the V 3d and Ti 3d orbital, respectively. Nb02 films are grown on Ill-oriented SrTiOa 
substrates as described in more detail elsewhere 23j. Both SrVOa and LaTiOa films are grown on 100-oriented 


SrTiOa substrates at a temperature of 600-800°G using co-deposition of matched metal fluxes in the presence of 
between 3 x 10“® to 2 x 10“® Torr of molecular oxygen with a total growth rate of ~0.4 nm/min. All films reported 
here are crystalline as-deposited, with pseudo-rutile structure for Nb 02 H and perovskite structure for SrVOa and 
LaTiOa, as determined by reflection high energy electron diffraction (RHEED). We systematically vary the oxygen 
pressure during growth to determine the conditions that would result in the ideal 0:Nb, 0:Ti, and 0:V ratios in the 
films. The transition metal to oxygen ratios are determined by the integrated intensities of the relevant XPS core level 
spectra (O Is for oxygen) and the appropriate atomic sensitivity factors, as well as verifying that the Sr:V and La:Ti 
ratios are very close to one. The atomic sensitivity factors used are empirical values as reported by Wagner et al. 
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2J,[23 and adjusted to give ideal oxygen to metal ratios for the compounds Nb 205 , V 2 O 5 , and undoped SrTiOa. In 


the following, we present our experimental results for the transition metal core level spectra for single phase, nominally 
d} materials, and for n-doped SrTiOs, as measured using in situ XPS. All materials are sufficiently conductive at 
room temperature such that there is negligible (< 0.1 V) sample voltage during the measurement. For each material, 
we show a core level spectrum for an under-oxidized, optimally oxidized, and over-oxidized sample for comparison. 
The detailed results for each material are presented in the following sections. In the Supplementary Material H , we 
provide RHEED data for stoichiometric SrVOa, LaTiOs, and Nb 02 to further demonstrate our sample quality. 


A. Nb02 

The Nb 3d core level in Nb 205 is located at a binding energy of 207.7 eV (Nb 3 d 5 / 2 ) and has a spin-orbit pair at 
2.7 eV higher binding energy (Nb 3 d 3 / 2 ). To model the Nb 3d multi-peak structure in Nb 02 , we assume that the 
spin-orbit pairs are of the same width and that their separation is the same as in Nb 2 05 . Two or three pairs of peaks 
(pseudo-Voigt line shape) are used as needed to fit the data. For the optimally oxidized case (0/Nb = 2.0) shown in 
Fig. [2] (a), we find two components. The first component has a binding energy of 206.5 eV with a width of 0.9 eV, 
while the second component has a binding energy of 207.5 eV and a larger width of 1.8 eV. If we assign the 207.5 eV 
feature to be the dP component, the dP component is 55% of the integrated intensity while the dP component is 45% 
of the integrated intensity. 

If we electron dope the system by removing oxygen to form an under-oxidized Nb 02 phase [Fig. [5] (b)] with 0/Nb 
= 1.9, we find that both the dP component at 207.3 eV and the d^ component at 206.0 eV decrease slightly in relative 
amount to 52% and 40% of the signal. A new component (d^) emerges at a binding energy of 204.5 eV with a relative 
amount of 8 %. On the other hand, if we add excess oxygen to the system and form over-oxidized Nb 02 [Fig. [2] (c)] 
with 0/Nb = 2 . 1 , the shape of the spectrum changes qualitatively. The dP component (at 207.6 eV) becomes sharper 
(width of 1.5 eV) and increases to 62%, while the d^ component at 205.9 eV (width of 1.1 eV) drops to 38%. 


B. SrVOs 

For SrVOa, we look at the V 2p core level. For comparison, in pure V 2 O 5 , the V 2 ^ 3/2 peak is located at a 
binding energy of 517.9 eV, with the 2 pi /2 spin-orbit pair located at 7.4 eV higher binding energy. The 2pi/2 peak 
is significantly broader than the 2p^/2 peak due to Coster-Kronig transitions. To model V 2p spectra, the widths of 
all 2 p 3/2 components are constrained to be the same and the widths of all 2 pii 2 peaks are also constrained to be the 
same. There is no restriction on the relative widths of the 2p^/2 and 2pi/2 peaks within each component, however. 
The 2 p 3/2 to 2pi/2 separation of each component is also fixed to be the same as that of V 2 O 5 . Three sets of spin-orbit 
pairs of peaks are used to fit all the SrV 03 data. Because the O Is core level is near the V 2p levels, O Is signals are 
also collected in the same measurement and included in the fitting. 

For the optimally oxidized case with 0/V = 3.0 [Fig. [3] (a)], the spectrum consists of three distinct components. The 
widths of the 2p^/2 peaks are 1.5 eV. The first peak has a binding energy of 517.9 eV {d^) with a relative concentration 
of 60%. The second peak (d^) has a binding energy of 516.2 eV with a relative concentration of 27%. The third peak 
{dP) has a binding energy of 514.5 eV with a relative concentration of 13%. Reducing the 0/V ratio to 2.7 [Fig.[3](b)] 
results in a significant decrease in the dP component at 518.1 eV to 36%. The dP component at 516.2 eV increases to 
36% while the dP component at 514.7 eV increases to 28%. On the other hand, slightly over-oxidizing the SrVOs to 
have an 0/V ratio of 3.1 [Fig. [3] (c)] alters the relative amounts of the three components to 65% for d°, 22% for d^, 
and 13% for dP. 
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C. LaTiOs 

For LaTiOa, we use the Ti 2p core level. The Ti 2pi/2 level is significantly wider than the 2p^/2 level due to 
Coster-Kronig transitions. We model the Ti 2p spectra using the same kind of constraints on widths and spin-orbit 
separation as in the V 2p modeling. For comparison, the Ti 2p^/2 level of stoichiometric SrTiOa (Ti"^+) is located at 
458.9 eV with a 2 ^ 3/2 to 2pi/2 separation of 5.6 eV. Two or three pairs of peaks are used to model the LaTiOa Ti 2p 
spectra as needed. For the optimally oxidized sample with 0/Ti = 3.0 [Fig. [4] (a)], there are two components. The 
first one (d°) is located at a binding energy of 458.4 eV with a relative concentration of 45%. The second component 
{d}) is located at a binding energy of 456.9 eV with a relative concentration of 55%. The widths of both 2 ^ 3/2 peaks 
is 1.7 eV. 

When LaTiOs is under-oxidized to yield an 0/Ti ratio of 2.8 [Fig.|4](b)], we see the emergence of a third component 
{(P) with a binding energy of 454.9 eV and a relative amount of 8%. The other two components are both slightly 
reduced in amount to 40% for dP and 52% for d}. For the slightly over-oxidized case, with 0/Ti = 3.1 [Fig. |4| (c)], 
we see a significant increase in the component to 69% with a slight shift in binding energy to 458.9 eV. The d^ 
component (at binding energy 457.0 eV) correspondingly decreases to 31%. 


D. n-doped SrTiOa 


In stoichiometric SrTiOs, the 2p^/2 peak shows a single feature about 1 eV wide with no shoulder [27|. The 2pi/2 


31 
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peak is significantly broader than the 2 ^ 3/2 peak due to Coster-Kronig transitions [2^. Fig. [^ shows the Ti 2 b XPS 
spectra for 15% La doped SrTi 03 (S ri_.T,La xTi 03 ) 
oxygen-deficient SrTi 03 (SrTiOs-x) 


19|,|29|, 10% Nb doped SrTi 03 (SrTii.^Nb^^Oa) [30|, and 


33| . In all these n-doped SrTiOs, a small shoulder located about 1.5 eV 


lower than the Ti"^+ peak emerges, and is typically interpreted as a Ti^+ (d^) peak. Two important features should 
be pointed out. First, the position and strength of Ti^+ peak are not sensitive to photoelectron emission angle (not 
shown), indicating that this signal is not a surface effect. Second, the position of Ti^+ peak is dopant-independent, 
indicating that this peak is very likely to be intrinsic to doped SrTiOs. We will show in the next section that these 
two observations are consistent with the final-state interpretation. 


E. Common features of transition metal oxide spectra 

We summarize this section by pointing out the key common features of the XPS spectra of these d^ transition 
metal oxides: the transition metal core level spectra of these materials all display at least two, and sometimes 
three distinct components (where a component refers to a pair of peaks related by spin-orbit coupling); a single 
component is never observed even in the optimally oxidized single phase films. These XPS peaks can be assigned 
as (Nb®’*', V®"*", Ti^+), d^ (Nb^"*", Ti^+), and d^ (Nb^’*', Ti^+) oxidation states. As a general trend, 

electron doping (via oxygen vacancies) increases the intensity of the d^ peak at the expense of the and d^ peaks, 
whereas hole doping (via oxygen excess) increases that of the peak and decreases the intensity of the d^ peak if 
present. Based on the initial-state effect, one might naively infer from the XPS results that the optimally oxidized 
samples contain significant amounts of regions of different oxidation states (such as Nb 205 which is nominally d°). 
However, this interpretation is not consistent with RHEED from the samples, which should clearly show the presence 
of incommensurate monoclinic/amorphous Nb 205 or pyrochlore La 2 Ti 207 /Sr 2 V 207 phases, if they are present in 
such large amounts. Quantitatively, if we assume the peak intensity of a particular component is proportional to 
the abundance of that particular oxidation state, this implies that roughly one half of the sample on average is 
in the highest oxidation state. For example, from the XPS of SrV 03 [Fig. [3], one expects 60% of the sample to 
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consist of pyrochlore Sr 2 V 207 which should be, but is not, reflected in the diffraction data, which still shows a 
single phase, epitaxial 100-oriented pervoskite film. It should also be noted that the oxygen to transition metal ratio 
has been carefully controlled during growth (as described above), spanning the range from under-oxidized to over¬ 
oxidized. Furthermore, we also note that growing at very low oxygen pressures that result in an oxygen to metal ratio 
significantly less than the ideal value still results in the presence of a peak that is associated with the dP oxidation 
state. The presence of a strong cf peak in stoichiometric SrVOa has been interpreted by Takizawa et al. 34 


as being due to excess oxygen (forming V®+) decorating the surface of SrVOa resulting in a ^2 x reconstruction 
pattern. As shown in the Supplementary Materials [2^, we also observe the surface reconstruction in RHEED. By 
comparing the XPS spectra before and after the Ar sputtering (which removes the surface atoms), we conclude that 
both the surface reconstruction (initial-state effect) and final-state effect contribute to the multi-peak structure in the 
case of SrVOa. In a vacuum-cleaved single crystal of SrVOa, only a weak dP feature is observable [Sg. The seemingly 
conflicting results from the XPS data and the single phase nature of the optimally oxidized films can be naturally 
reconciled if the occurrence of the multi-peak structure in the XPS spectra is intrinsic to these d^ materials (i.e. the 
spatially uniform d^ system by itself displays multiple peaks in XPS). In the next section we argue that it is indeed 
the case once the final-state effect is considered, and provide a simple model to illustrate this point. 


IV. MODEL AND THEORETICAL ANALYSIS 


A. Model and parameters 


To explain the observed multi-peak structure in the XPS spectra, we propose a cluster-bath model which resembles 

|. It contains three parts [see Fig. |7Ka)]: 


that proposed in Ref. 


37 


Htot — ^cluster T dIbath did — bath • 


(3) 


Hduster describes a TM-Og (TM can be Ti, V, or Nb) cluster that includes at least five TM d and three O 2p 
orbitals (for each of six oxygen atoms). We point out that our model does not qualitatively distinguish between the 
3d and 4d orbitals, or 3d orbitals of different chemical elements: they only correspond to different parameters in the 
model. When taking the cubic symmetry into account, only ten of twenty-three total orbitals couple to one another 


39j. Using F to label the orbital symmetry (three t 2 g xy, yz, zx and two eg 'iz'^ — r^, x'^ — y^ orbitals), the cluster 


Hamiltonian 


401 is 


Hduster — ^ ^ Cp(F)7rp^r,o 
r,a 


ed(r)nd,r.cr + U(F)[djf^^pr,o- + id.c.]| 


-f 


U 


2 S 


'bld^r.a'bld^r' Uic(l 'btcoreP ^ ^ 'bldp^a “t" Cc 


(4) 


Here nd,r,a- = dp^dr.o-, ^^p.r.o- = Pp crPr.o- are respectively the TM d and O 2p number operators for the orbital 
labeled by (r,cr) (cr labels the spin). ed(r) and ep(r) are energies of TM d and O 2p orbitals, and U(r) describes 
their hybridizations. U is the energy cost when the 3d occupation of the TM atom is more than one. ricore is the 
number operator of the core level and Cc the core-level energy (approximately —459.0 eV for Ti 2p, —518.0 eV for V 
2p, —208.0 eV for Nb 3d). The term with Udc approximates how valence electrons respond to the core hole: in the 
presence of a core hole (i.e. (ricore) = 0), all TM d levels are shifted down by Udc to screen the core hole. By fitting 
to published experimental data from XPS of SrTiOg, ellipsometry, and angle-resolved photoemission spectroscopy 


(ARPES) [ia,l32|,l4ll-l43i, we take edicg) = 2.0 eV, ed(t 2 g) = 0 eV, ep(r) = -3.0 eV, V(eg) = 2.5eV, V(t 2 g) = -1.3 
eV, U = 6.0 eV, and Udc = 8-0 eV This problem can be solved exactly by the technique introduced by Gunnarsson 
and Schonhammer 


12 


45l| . and the details are provided in Appendix B. 
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For the metallic phase, the occupation of each local orbital fluctuates. To capture this effect, we further introduce a 
set of bath orbitals, which simulate the role of TMO conduction bands, coupling to each d orbital. For Hamiltonians 
involving the bath: 


Hbath Hcl — bath — ebl-^^bercrde + J [Vie,T)dl^bercr + H.c.jde, 


2V^ 


7r|H(e,F)p = 


(5) 


Here b^ro- denotes the bath orbitals of energy e, orbital symmetry F and spin cr. Inclusion of the bath introduces 
charge fluctuation in the cluster (via exchange of particles with the bath) that is used to model the fluctuation in the 
occupation of local orbital in the metallic phase [l2| (see Appendix A for a simple explanation). We use ep = 2.0 eV, 
B = 2.0 eV (so the bath levels range from 0 to 4.0 eV, roughly the SrTiOp conduction bandwidth) to approximate the 
SrTiOp conduction bands, and take V = 0.3 eV which is approximately the effective hopping between two adjacent 
Ti 3d orbitals 4l|, l46|. It turns out that the exact value of V plays a relatively minor role in the XPS spectrum (see 


Appendix B). In the calculation, we introduce the chemical potential /r to specify the number of total electrons (filling) 
in the whole system (bath and cluster): all bath levels below fi are filled. Qualitatively larger fi corresponds to larger 
average d occupation in the bulk material. To extract the essential feature of these d^ materials, we only vary /r but 
keep all other parameters fixed. In other words, the valence levels of Ti 3d, V 3d, and Nb 4d are not distinguished in 
our simulation. The XPS spectrum is calculated using Eq. and the details are given in the Appendix B. 


B. Results from an isolated cluster 


Before discussing our results using the total Hamiltonian Eq. ©, we first present the results from the isolated 
cluster (zero impurity-bath coupling). In particular we shall identify the origin of each peak. Eig. [IKc) shows the 
XPS spectrum of an isolated cluster - ten electrons are filled to mimic the nominally d^ system. There are three 
pronounced peaks, labeled as \L), |M), and \U) referring to their relative lower, middle, and upper bindin g en ergies. 
These features can be understood by considering the following three states |d°L°), |d^L^), and |d^L^) [s, 15|. Here 
|d°L°) represents the “reference” state where all O 2p orbitals are filled, and represents the state of i particle- 

hole (p-h) pairs with respect to [see Pig.HKa) for illustration]. Without the core hole, the ground state \GS) is a 

linear combination of these three states. States with larger number of p-h pairs are significantly less important due to 
the on-site energy U. In the presence of a core hole (we use |d*T*c) to denote states in the presence of a core hole), the 
relative energies of these three states change, and the resulting core-hole eigenstates (including the d-p hybridization) 
are labeled as |L), |M), \U). These three lowest eigenstates account for the three pronounced peaks in the computed 
spectrum. From Eq. (ED, the peak strength is given by |(A|c|G'5')p for X = L,M,U. We emphasize that, due to 
the strong d-p hybridization, all core-hole eigenstates |L), |M), \U) have significant |d*T*c) (z = 0,1,2) components. 
Comparing with the experimentally observed XPS SrTiOa spectrum [Eig. [5], we note that: (i) the strongest peak 
\L) is conventionally assigned as the Ti^+ (d^) 2p^/2 peak; (ii) the weak peak \M) is buried under the 2pi/2 peak 
caused by the spin-orbit coupling of the core electron, and is not observed; (iii) the calculated \U) peak corresponds 
to the charge transfer satellite feature at a binding energy of approximately 471.0 eV 15, ^ and appears to be much 


sharper than that in the experiment, because we neglect the coupling between valence electrons and the core spin 
that provides additional decay channels for states of higher binding energies 15, 4^ . In the following discussion we 


only focus on the strongest and lowest peak, which is the one used to determine the different oxidation states. 















C. Results including bath 


Inclusion of the bath introduces charge fluctuation in the cluster, as the cluster can now exchange particles with 
the bath orbitals (see Appendix A). More specifically, instead of the fixed number of electrons in the cluster, the total 
ground state wave function has a general form 

^ ^ '^)cl ^ 1^6 ‘^)bath ^ |l)core- (6) 

i=0 

Here \nci + i)ci ® \nb — i)bath ® |l)cor-e represents a state which has rici + i particles in the cluster, rib — i particles in 
the bath and a filled core level. Using the same notation as Eq. ®, the isolated cluster calculation presented in the 
previous subsection only has \nci = 10 )ci ® \nb)bath ® |l)cor-e, with \nci = 10 )c/ including all possible (f =0 to 10 

in principle) components. When exchanging particles with the bath, the states such as \nci = ll)ci * = 0 

to 9), \nci = 12)cj f = 0 to 8 ) also contribute to the \GS). Similarly, in the presence of a core hole, the nth 

eigenstate with energy Ecore,niN — 1 ) has the general form 

\n{N - 1)) = ^ \nci + i)ci 0 \nb - i)hath 0 |0)core- (7) 

Applying Eq. ([ 6 ]) and Eq. 0 to Eq. 0 , the XPS spectrum displays peaks at = Eas{N) — Ecore,n{N — 1) with 
weight I Erom this general analysis, we see that including the charge fluctuation naturally leads to 

multiple XPS peaks, which correspond to different particle number in the cluster. 

In Fig. [7Kc) we show the calculated XPS spectra for /r = 0.2,1.0,1.5, 2.0, and 2.5 eV. Starting from the highest 
chemical potential, the /i = 2.5 eV XPS spectrum shows three distinct peaks. By analyzing the wave functions, they 
correspond to |ncz = 10)cz, |ncz = ll)cZ and |ncz = 12)cz in Eq. ©, and are therefore labeled as d}, (P respectively. 
Using the notation within the isolated cluster, the fi°, d} and peaks come from states of \L) (s |ncz = 10)cz), 

(G |ncz = ll)cz), (g |ncz = 12)cz) respectively, as shown in Fig. [7Kb). Decreasing p, reduces the intensities of 

and d^ peak but increases that of <P. This is because lowering the chemical potential decreases the probability of 
adding electrons to the cluster from the bath, resulting in a smaller jricz = ll)cZ and jricz = 12)cz components in \GS) 
and consequently weaker d^ and peak intensities. However, we stress that once the cluster and bath can exchange 
particles, a single d^ XPS peak is never obtained in our calculation; the d^ peak is always present. 


D. Discussion 


1. Comments on experiments 


We now discuss several experiments based on the calculation. The main conclusion from our model calculation 
is that, once charge fluctuation is taken into account, the nominally d} or n-doped dP transition metal oxides are 
expected to display multiple peaks in their XPS spectra, even in the absence of other oxidation states. In other 
words, our theory implies that a multi-peak structure in XPS is general for these materials if charge fluctuation 
cannot be neglected. 

We first discuss three observations in n-doped SrTiOa samples based on the general consequences of the final-state 
interpretation. First, the Ti^’*' peak position is dopant independent and is an intrinsic property of the Ti atom, or 
more precisely the TiOg cluster . Indeed, in lightly n-doped SrTiOa, the Ti^’*' peaks all appear in the same position 
relative to the TU+ peak 


17 


19(1 (Fig. [5]). Special attention is paid to the Nb-doped SrTiOa (or Nb-doped Ti02 [49|), 


where even in the ionic limit, there can only be Nb^’*' ions (i.e. Nb keeps one Ad electron), but not Ti^+ ions. Within 
our interpretation, the Nb gives its 4fi electron to the conduction band, resulting in a metallic state and nominally 
Nbl®-^)’*' and Ti(^““)+ ions (instead of Nb'*’*' and TU+), with Ti(^“^)+ ions providing the XPS TU+ signal. Second, 
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the Ti^+/(Ti"^++Ti^+) ratio is routinely used to estimate the dopant concentration, and gives very reasonable values, 
which are consistent with other experiments such as Hall measurements and Rutherford backscattering for low to 
moderate doping 17 J3- According to our theory, this is possible because the Ti'^+ is the highest oxidation state 
and contains only one main peak. The Ti^+ signal therefore appears as an extra, distinct side peak when reducing 
the average Ti oxidation state via doping. For a nominally d} system (that will be discussed shortly), multi-oxidation 
peaks exist intrinsically in the first place, and doping does not introduce a new peak. Also, we expect that using 
the Ti^+/(Ti"‘+-|-Ti^+) ratio always slightly underestimates the dopant concentration as the nominally pure Ti^+ 
material already has significant Ti"‘+ signal. This is consistent with the results in Ref. [l^. Finally, one cannot really 
distinguish the initial-state and final-state effect based solely on the XPS spectrum. Both spatially localized Ti^+ 
ions or a uniformly distributed Ti(^““^+ can account for the XPS Ti^+ peaks. The key difference between these two 
scenarios is that the former implies the presence of an in-gap state, whereas the latter does not. To differentiate 
between them, one should probe the valence states to see if there is an in-gap signal. In oxygen-deficient SrTiOa, 
an in-gap signal is observed in ARPES 31, 32, In this case the XPS Ti^+ peak can be due to the presence of 


localized Ti^+ ions. We note that in the literature, an oxygen vacancy is suggested to be a single donor l52l| . 


which would result in nominally localized Ti^ ®+ ions (we favor this view). Within the final-state effect, Ti^-®+ ions 
also lead to a separate XPS Ti^+ peak. It is worth no ting that in the LaAlOa/SrTiOa interface, the oxygen vacancies 


55, 


are responsible for the majorityof charge carrier 53l. l54|. However, the x-ray absorption spectrum does not indicate 
the existence of Ti^+ ions 


56|. 


For the nominally TMO, all the optimally oxidized d^ samples we have grown (as well as vacuum-cleaved single 
crystal Ti 203 ^), demonstrate the XPS spectra showing multiple components. By viewing the multi-component 
structure as being caused by the final-state effect, the existence of these multiple components does not require the 
presence of different oxidation states in the sample. Even though XPS data show multiple components, the systematic 
way in which the oxygen content is controlled in the growth experiments, in combination with the single phase RHEED 
patterns observed, precludes the existence of different oxidation environments in the optimally oxidized samples. The 
final-state interpretation reconciles the seeming conflict between XPS data, the single phase pattern in RHEED 
measurements, as well as the careful, systematic way in which the oxygen content is controlled in these growth 
experiments, which precludes the existence of different oxidation environments. Moreover, our calculation shows the 
same doping dependence of the relative peak intensities: increasing the electron doping decreases the d^ peak intensity 
and causes an increase in the intensity of the peak. This qualitative agreement between theory and experiment 
leads us to believe that the multi-peak structure in the single phase d^ transition metal oxides actually originates from 
the final-state effect and is intrinsic. Certainly, as mentioned previously, one cannot rule out the initial-state effect, 
and ions of higher oxidation states (F®^ for example) may exist at or near the surface of the sample. As observed 
in some vanadates 34-^, these ions also result in d^ signal. However, we notice that even if these ions do exist, 


the d^ signals appear to be too strong (d° and d^ peaks are of comparable strength) to be interpreted as being solely 
from them. In fact, we believe in SrVOa, the surface reconstruction (initial-state effect) and final-state effect both 
contribute to the observed d° peak (see the Supplementary Materials [26|). 


2. Limitations of the theory 

There are two uncertainties in our model which make a more quantitative analysis difficult. Eirst it is not easy to 
map the chemical potential /i to the average d occupancy in the bulk material. Second, the energy distribution of bath 
orbitals and the cluster-bath coupling are also hard to determine. However, the multi-peak structure is insensitive 
to these uncertainties (see Appendix B). Namely, as long as there are particle exchanges between the cluster and the 
bath, there are multiple peaks in the XPS spectrum. For this reason we believe the conclusions drawn from our model 
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are qualitatively correct. 


3. Charge fluctuation 


We now discuss the origin and the importance of charge fluctuation. The charge fluctuation cannot be neglected in 
the metallic state, where particle exchange with the Fermi sea causes fluctuation in the occupation of local orbitals 
[l^ . Accordingly, charge fluctuations in doped or metallic samples should not be neglected, and multiple XPS peaks in 
these samples are expected (and indeed observed) [^. For undoped, nominally insulating d} materials, the criterion 
of being metallic is not always satisfied at low temperature. Uncorrelated d} materials are expected to be band 
metals. The samples we have studied, Nb 02 and LaTiOfj, are both metallic at high temperature and undergo a 
metal-to-insulator transition at 1080 K (of Peierls type) [57j and 125 K (of Mott type) [58| respectively; SrVOa is 
intrinsically metallic 0 . Note that SrVOa already shows the d^ peak in the optimally oxidized sample [Fig. [3] (a)], 
indicating its relatively strong charge fluctuation due to its metallic nature. 

Specific to our experimental conditions, all n-doped SrTiOa are metallic at room temperature. LaTiOa is already 
metallic at room temperature, which easily allows for charge fluctuation. For Nb 02 , the sample is still nominally 
insulating at room temperature, but its relatively small band gap of ^1.0 eV 0| likely results in non-negligible 
concentration of electrons in the conduction band at room temperature. The fact that no sample charging is observed 
during XPS measurements indicates that there is sufficient conductivity in the samples at room temperature (sufficient 
thermally excited carriers in the conduction band) to allow for charge fluctuation to occur. Therefore, although the 
charge fluctuation in the undoped, nominally insulating dd materials can be weaker compared to the doped samples, 
we believe it is still non-negligible. 


E. Relative importance of the initial-state and final-state effects 


We would like to conclude our theoretical analysis by addressing the relative importance of the initial-state and 
final-state effects. From Eq. 0, we see that the valence screening is described by the parameter Udc-, which is the 
strength of the core-hole-induced attractive potential. If Udc = 0, then valence band electrons do not feel the existence 
of the core hole, and thus no final-state effect is involved. With this observation, we propose that the dimensionless 
parameter ^ = Udc/W^ with W the typical energy scale of the valence bandwidth, can be used to characterize the 
relative importance between initial-state and final-state effects: large ^ favors the final-state effect; small ^ favors 
the initial-state effect. As the bandwidth is proportional to the electron hopping f, we can roughly regard IjUdc 
as the time scale to create a core hole, and 1/t as the time scale for conduction electron to move to screen the 
core hole. Therefore the inverse of ^ (1/C) essentially describes how efficient (fast) the conduction electrons screen 
the core hole. By fixing the value of Udc (about 10 eV [^), materials of large/small valence bandwidth favor the 
initial-state/final-state effect. 

With this picture, we comment on the established interpretations of XPS spectra. For covalent materials such as 
carbon and silicon, the initial-state appears to be dominant and the multi-peak structure is used appropriately to 
signal the existence of different oxidation phases 0,01 ■ Consistent with our ar gum ent, the diamond structure of C 
and Si indeed have relatively large valence bandwidths of approximately 20 eV 60, 0 and 12 eV 0|, respectively. 


which favors the initial-state effect. For materials with valence electrons in localized orbitals (rare earth 4/) such as 




lanthanum and cerium [ill,[ij,l22| , it is the final-state effect which dominates. For these materials the XPS multi-peak 
structure is not attributed to the oxidation states, but can be used to determine material-specific model parameters 
by comparing to a model calculation 0 . A typical bandwidth of /-orbitals is about 4 eV 21, 63, 0 , which favors the 
final-state effect. In terms of valence bandwidth, the early transition metal oxides are in between the two classes of 
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materials (about 6 to 8 eV 4l|, [4^) . As the experimental results for carefully grown samples from different probing 


techniques fit the final-state effect better (see also Refs. 37|,l3^), we believe the final-state effect is also the dominant 


one in the transition metal oxides. Taking Udc to be 10 eV, we summarize the origin of the multi-peak structure in 
XPS for the materials mentioned above in Table ID 


Material 

valence bandwidth {W) ^ = 

UdclW origin of multiple peaks 

diamond carbon 

20 eV 

0.5 

initial-state 

diamond silicon 

12 eV 

0.83 

initial-state 

SrTiOs 

6 eV 

1.67 

final-state 

CeNiz 

4 eV 

2.5 

final-state 


TABLE I: The origin of the multi-peak structure in XPS for various materials. The value SrTiOs is similar to the materials 
studied in this paper. Udc is taken to be 10 eV. 


V. CONCLUSIONS 

We investigate the origin of the observed XPS multi-peak structure of single phase nominally transition metal 
oxides including Nb02, SrVOa, LaTiOa, and lightly n-doped SrTiOs. Experimentally, we find that the XPS spectra 
(specifically the photoelectrons from Nb 3d, V 2p, Ti 2p core levels) of these materials all display at least two, and 
sometimes three pairs of peaks, which can be consistently assigned as d^, and d^ oxidation states. For lightly 
n-doped SrTiOs, a weak d^ shoulder, whose energy position is independent of the dopants, appears with respect to 
the main peak. For nominally d^ transition metal oxides, electron doping increases the intensity of the peak 
but decreases that of the d° peak, whereas hole doping reverses this trend. A single d^ peak is never observed, 
even in single phase samples. In particular, the d^ peak always exists even in the electron doped samples where 
stoichiometric analysis shows strong oxygen-deficiency and diffraction shows no secondary phases, strongly indicating 
that the multi-peak structure is intrinsic to these materials. Theoretically, we construct and solve a cluster-bath model, 
and explicitly demonstrate that the final-state effect (i.e. the valence response to the created core hole) naturally 
leads to the multiple peaks in the XPS spectrum even in a spatially uniform system. Moreover, the relative peak 
strength as a function of doping is qualitatively consistent with the experimental observation. The combination of 
experimental and theoretical analysis leads us to conclude that the multi-peak structure in the nominally d^ transition 
metal oxides is intrinsic, and does not necessarily imply the existence of spatially isolated (or clustered) d*^ and d^ 
ions in a sample. Using the same analysis, we argue that the ratio between the local screening potential and the 
valence bandwidth is the key dimensionless parameter that determines the relative importance between initial-state 
and final-state effects. To establish the existence of different oxidation phases in a sample, further spatially-resolved 
probing techniques involving the valence electrons are needed. For this reason, investigating the final-state effect in 
x-ray absorption spectroscopy can be very helpful. 
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Appendix A: Charge fluctuation in the metallic phase 

In Section IV we propose a model [Eq. ([5])] and argue the importance of charge fluctuation in the metallic phase. 
Here we use a very simple model to illustrate this effect. We consider a three-site tight-binding model containing two 
electrons: 


H 30 -TB = + C2,<tC3.<t + + h.C.), 

<7 


(Al) 


with a the local orbital basis. The two-particle ground state is |(/>o) = dj ^\vac) with = (c| -I-C2 -I-C3 g.)/\/3. 

Note that dj describes a Bloch orbital which is spatially extended. When expressing the ground state using the 
local orbital basis, we have 


\^o) = ^c\^^c\^^\vac) + ^ 


3 1 


.t 


(4,4. + C' 


) + (< 


'2,t ■ 


4,4'' 


C4 ..)c 


\vac) + 3 ( 4 , 4 ' 


l)(' 


3 , 4 '''' 2,4 


"3,4'' 


(A2) 


The terms are grouped according to the occupation on the first site. In the local basis, we see that the many-body 
ground state (only two-body in this case) contains all doubly occupied, singly occupied, and unoccupied components; 
the situation which is typically referred to as the charge fluctuation. Note that this fluctuation has nothing to do with 
temperature, but originates solely from the many-body wave function. The inclusion of the bath degree of freedom is 
to take this charge fluctuation of the local occupation into account. 


Appendix B: Details of computing core-level spectra 


1. Basic formula 


The XPS spectrum is computed using 



p(cc) = ^ |(n{iV - l)|c|G5)p X d(a; - [Egs{N) - E,ore,n{N - 1)]) 

n 


= 4(G5|ct [w - {EasiN) - Htot) - id]”' c|G^). 

TT 

In the calculation, we first solve the ground state |G5'), and then follow the Lanczos procedure 
the spectrum. 


(Bl) 


65 


to compute 


2. Cluster model 


Here we compute p{uj) for the cluster model specified in Eq.(2) in the main text. This problem can be exactly 
solved b y di agonalizing only a 35 x 35 matrix. The key numerical step, first realized by Gunnarsson and Schdnhammer 
in Ref. [1^, is that for all degenerate determinantal states, only one of their combinations contributes to the exact 
ground state. In Table HIl we list all 35 states. The reference determinantal state, |0), is defined by occupying all O 
2p levels, and the other 34 states are labeled by particle-hole (p-h) pairs in t 2 g and eg sectors. 

To explicitly write down these states, we define the p-h operators as Pi = d\pi for i belongs to one of six t 2 g orbitals 
(including spins), and Pj = d^^pj for j belongs to one of four eg orbitals. The state labeled as {nt;nt 2 g,neg) in Table 
His 


(n^, 72^2(7, ^ 



1 

w 

1 _ 





.G> 


| 0 ). 


(B2) 
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p-h pairs (notation) 

label 

{nt', nt 2 g , Ueg) [u] 

0 (|d°L°)) 

|0) 

(0; 0, 0) [1] 

1 (IdiLi)) 

|1) to |2) 

(1; 1, 0), (1; 0, 1) [2] 

2 {\d^!?)) 

|3) to |5) 

(2; 2, 0), (2; 1, 1), (2; 0, 2) [3] 

3 (|d^T"» 

|6) to |9) 

(3; 3, 0), (3; 2, 1), (3; 1, 2), (3; 0, 3) [4] 

4 (Id^L^)) 

|10) to |14) 

(4; 4, 0), (4; 3, 1), (4; 2, 2), (4; 1, 3), (4; 0, 4) [5] 

5 (|d®L®)) 

|15) to |19) 

(5; 5, 0), (5; 4, 1), (5; 3, 2), (5; 2, 3), (5; 1, 4) [5] 

6 (|d«L®)) 

|20) to |24) 

(6; 6, 0), (6; 5, 1), (6; 4, 2), (6; 3, 3), (6; 2, 4) [5] 

7 (|d^L^)) 

|25) to |28) 

(7; 6, 1), (7; 5, 2), (7; 4, 3), (7; 3, 4) [4] 

8 (|d«L«)) 

|29) to |31) 

(8; 6, 2), (8; 5, 3), (8; 4, 4) [3] 

9 (|d^L®)) 

|32) to |33) 

(9; 6, 3), (9; 5, 4) [2] 

10 (M“L^°)) 

|34) 

(10; 6, 4) [1] 


TABLE II: States in the cluster with ten electrons, nt, nt 2 g, n^g are the number of total p-h pairs, t 2 g p-h pairs, and Cg 
p-h pairs respectively, [n] represents the number of states for the number of p-h pairs. The diagonal energy is given by 
Unt{nt - l)/2 - nttp + neged{eg) -\- nt 2 gtd{t 2 g)- 


Here {z} ({j}) represents all combinations of creating nt 2 g {rieg) p-h pairs out of the reference state, and C:^ = 
(jv-^)!m! ■ that each individual determinantal state in the summation has the same energy, and it is Gunnarsson 
and Schonhammer’s invaluable observation that only the sum of them contribute to the exact ground state, and all 
other combinations can be rigorously neglected. The coupling between these states is nonzero only when |Aneg| -|- 
= 1, and can be computed straightforwardly. To check this formalism, we also computed the ground state 
and XPS spectrum by diagonalizing the original 63504 x 63504 matrix (the dimension of filling 10 electrons in 20 
orbitals with = 0), which gives identical results to that obtained by keeping only 35 states. When including the 
bath degrees of freedom, it is not possible to include all states. In Fig. [8] we show that keeping states up to two p-h 
pairs already results in a very reasonable profile. Keeping states up to four e-h pairs almost reproduces the exact 
spectrum. 


3. Cluster coupling to bath 


Now we solve the problem including both Eq. (j4]) and Eq. ([5]) in the main text. There are six degenerate t 2 g and 
four degenerate Cg d-p orbital pairs in the cluster, and each 3d orbital couples to its own bath. In the calculation, 
it is convenient to simply treat the O 2p orbital as one of the baths, whose energy and coupling to the Ti 3d are 
respectively ep(r) and K(r). In other words, we should notationally identify = Pr, and V{€p,T) = K(r). The 
reference state |0) is chosen as 


| 0 )^ 




\vac), 


(B3) 


i.e. all bath levels below the chemical potential p, are filled. 

Similar to the previous subsection, we define the p-h operators for the t 2 g and eg sectors: Pi{e) = d\b^^i Pi{E, e) = 
6 ^ -b^^i for i G t 2 g; Pi{e) = d\b^^i Pi{E,e) = 6^ -bf^^i for i G eg. The E and e are for bath states which are higher and 
lower than the chemical potential. We have tested energy spacing by discretizing the bath continuum into 32 to 100 
intervals, and they result in essentially identical spectra. We keep the states up to 2 p-h pairs, as tested in Ref. [l^ . 
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They are (in addition to the reference state) 

\ed)t2g 

I ^d) eg 
\f^E)t2g 
\iE)eg 
led, e'd)t2g 
\ed, ed)t2g 
\ed, ed)eg 
led, ed)eg 

|ed, e d)jriix 


1—1 

^^P.(£;,e)| 0 ), 

^ 2=1 

^^P.(£;,e)| 0 ), 

^ 2=1 


5]P,(e)P,,(e)|0), 


x/6(6-1)/2J^ 




6 4 




(B4) 


The coupling between two states is non-zero only if the number of p-h pairs differs by one. The XPS spectra are 
computed within these states. Finally, in Fig. El we show the computed XPS spectra for ^ = 1, and V =0.1, 0.2, 0.3 
eV (defined in Eq. ([5]) in the main text), with the Ti^+ peak appearing for all of them. We emphasize again that 
the main role of the bath coupling is to introduce charge fluctuation within the cluster, and the form of the coupling 
plays a relative minor role in the spectrum. 


Supplementary Materials 
RHEED data 

We provide in Fig. [10] the RHEED patterns for optimally oxidized epitaxial films of nominally d^ transition metal 
oxides: (a) LaTiOa; (b) SrVOa; (c) Nb02. For LaTiOa and SrVOa the RHEED pattern shows four-fold symmetry 
with a weak 2x reconstruction for LaTiOa and no reconstruction for SrVOa. For Nb02, the RHEED pattern shows 
six-fold symmetry due to the existence of three symmetry-related rotational domains. 


Surface effect: SrVOa XPS spectra before and after Ar sputtering 


Here we provide the SrVOa XPS spectra with and without Ar sputtering, to show that both the surface reconstruction 
(initial-state effect) and final-state effect contribute to the multi-peak structure. In situ XPS spectra of SrVOa thin 
films before and after sputtering SrVOa is thin films exhibit a weak y/2 x \f2 surface reconstruction in low energy 
electron diffraction 3^. This surface shows two components in the O Is spectrum and at least two components in the 


V 2p spectrum. This is consistent with an o xyg en-rich surface layer with species present, even when the growth 


and XPS measurement are done in situ 


34 l35 |. Angle-dependent XPS indicates that this oxygen-rich layer is 2-4 A 
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thick |34| . We performed in situ XPS measurements of a SrVOa thin film before and after argon ion sputtering of the 
surface to determine the effect of removal of the oxygen-rich surface layer on the O Is and V 2p spectra of SrVOa. 
A stoichiometric SrVOa film is grown by molecular beam epitaxy on a SrTiO3(100) substrate at a temperature of 
700°C to a thickness of 150 A. At this thickness, none of the substrate photoelectrons is able to escape the surface, 
allowing for a clean determination of the oxidation states and stoichiometry of the film. XPS analysis is performed in 
a separate analysis chamber that is connected to the growth chamber by a vacuum transfer line with a base pressure 
of 1 X 10“® Torr. The analysis chamber has a base pressure of 3 x 10“^° Torr. XPS spectra are taken before and 
immediately after sputtering. The sample is never exposed to the ambient at any time. The surface sputtering is 
performed using argon ions generated by a differentially pumped Hiden IG 20 ion gun. The Ar background pressure 
is 5 X 10“® Torr using an electron emission current of 20 mA. The Ar ions are accelerated to an energy of 1 keV and 
rastered on 3 mm x 3 mm area of the sample resulting in a sample current of 120 nA. The ion sputtering is performed 
for a total time of 5 min. Fig. ITlT al shows the O Is/V 2p spectrum before sputtering. 

The spectrum before sputtering can be deconvoluted into two oxygen components (531.0 and 529.7 eV) and three 
vanadium components. The oxygen feature at 531.0 eV corresponds to strong surface oxygen feature consistent with 
the presence of an oxygen-rich surface layer. The oxygen feature at 529.7 eV corresponds to bulk oxygen in SrVOa. 
The vanadium components have binding energies of 517.9, 516.2 and 514.5 eV corresponding to 3d° (V®+), 3d^ 
and 3(P (V^~*~) components. The strong 3d° feature is consistent with a surface layer as described by Takizawa 
et al. 3^. Note that the oxygen-rich surface reconstruction is present even when the sample has not been exposed 


to excess oxygen. The O Is/V 2p spectrum after sputtering is shown in Fig. [TTT al. The same set of two oxygen 
components and three vanadium components are still present after removal of the surface atoms in the film. However, 
the surface oxygen component is reduced by a factor of 3/4 showing that most of the surface oxygen has been removed. 
The 3dP vanadium component is at the same time reduced by more than half (with a slight shift in binding energy 
to 518.1 eV). This shows that even after removal of the surface layer, the 3(f component is still present although 
significantly reduced in magnitude. The after sputtering spectrum is now more like the bulk-sensitive SrVOa spectrum 
reported by Eguchi et al. using hard x-ray photoemission [^ . 

Because of the occurrence of an intrinsic surface reconstruction resulting in an oxygen-rich surface layer, which 
results in the presence of a 3fi° feature in the V 2p spectrum, it is not possible to ascribe the observed 3fi° feature as 
arising solely from the final state effect described in the main text. However, removal of the surface layer still shows 
the presence of a significant although strongly reduced 3dP component. While most of the 3cP signal in the as-grown 
SrVOa film is due to the surface reconstruction, there is a residual 3cP component in SrVOa present, which is what 
we ascribe to being due to the final state effect. 
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FIG. 1: (Golor Online) Illustration of XPS measurements, (a) (Left) The overall initial state contains an incoming photon of 
energy hv, and the system is in the ground state with the core level filled. (Right) The overall final state contains an outgoing 
electron of kinetic energy Ekin, and the system is in some excited state with the core level empty, (b) For the initial-state 
effect, the valence electrons are not affected by the created core hole, and therefore the binding energy is determined by the 
core-level energy Cc. (c) For the final-state effect, the valence electrons do respond to the created core hole. This can lead to 
multiple peaks around Cc in the XPS spectrum. 
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FIG. 2: (Color Online) Nb 3d XPS spectrum of Nb02: (a) optimally oxidized; (b) under-oxidized; (c) over-oxidized. The 
arrows indicate the dP, and d^ peaks. The S' jdP peak is most pronounced in over-oxidized/under-oxidized samples. 
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FIG. 3: (Color Online) V 2p XPS spectrum of SrVOs: (a) optimally oxidized; (b) under-oxidized; (c) over-oxidized. The 
arrows indicate the d°, d^ and d^ peaks. The d°/d^ peak is most pronounced in over-oxidized/under-oxidized samples. We note 
that in SrVOa, the optimally oxidized sample already displays a dP peak. 
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FIG. 4: (Color Online) Ti 2p XPS spectrum of LaTiOa: (a) optimally oxidized; (b) under-oxidized; (c) over-oxidized. The 
arrows indicate the dP, and d^ peaks. The S' jdP peak is most pronounced in over-oxidized/under-oxidized samples. 






FIG. 5: (Color Online) Experimental Ti 2p XPS spectra for stoichiometric, 15% La-doped, 10% Nb-doped, and oxygen- 
deficient SrTiOa. A shoulder, labeled as Ti®^, at about 1.5 eV below the main Ti'^'*’ peak appears for all n-doped sam ples . 
The dots are experimental data and the solid curves are from Gaussian fitting. The experimental details are given in Ref. [lOl ]. 
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FIG. 6: (Color Online) (a) The orbitals kept in the cluster. For the state \d°L^), all 2p levels are filled. For the state 
one of the 3d levels is filled, resulting in a 2p hole of the same orbital symmetry. For the state |d^L^), two of the 3d levels are 
filled, (b) A simple energy spectrum, keeping only the \cf LP), |d^L^), |d^L^) states. The double-arrow indicates the non-zero 
coupling between these states. The ground state \GS) is a linear combination of these three states. In the presence of a core 
hole (we add c to indicate its presence), the relative energy levels change due to the screening, but their hybridizations stay the 
same. The resulting eigenstates, labeled as |L), |M) and \U), correspond to the peak-structure in the XPS spectrum, (c) The 
computed XPS spectrum in arbitrary units, computed with a broadening of 0.2 eV. Three main peaks can be understood as 
the excitations in the presence of a core hole. The peak \L) is conventionally assigned as the d° peak (in the sense of |d°L°)), 
although it contains significant contributions from |d^L^) and |d^L^). 
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FIG. 7: (Color Online) (a) Schematic illustration of the cluster impurity model. Each d-orbital couples to a set of uncorrelated 
bath orbitals, (b) Cluster energy spectrum with a core hole. Due to the bath coupling, the number of electrons within the 
cluster is not a constant. In particular, the state of |d^L°c), |d^L°c) accounts for the d^, d^ peaks at lower binding energy, (c) 
The computed XPS spectra (arbitrary units) for p — 0.2,1.0, 1.5, 2.0, 2.5 eV. The spectrum of p = 0.0 is given as the shaded 
region. Increasing the chemical potential increases the intensity of d^ and d^ peaks and reduces the d° peak. For these chemical 
potentials both dP and peaks always exist. At p = 2.0 eV, the d^ peak begins to emerge as a shoulder at the lower binding 
energy. A Lorentz broadening of 0.2 eV is used. 
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FIG. 8: (Color Online) The p-h number dependence of the XPS spectra (arbitrary units) for the cluster model. The shaded 
area represents the exact spectrum. We see that including states up to two p-h pairs already gives a very reasonable profile. 
Keeping states up to four p-h pairs e 



FIG. 9: (Color Online) The computed XPS spectra for /r = 1.0 eV and V=0.1, 0.2, 0.3 eV. The Ti®^ peak emerges for all 
these V values. A Lorentz broadening of 0.2 eV is used. 



FIG. 10: (Golor Online) RHEED patterns for optimally oxidized epitaxial films of nominally transition metal oxides, (a) 
LaTiOa film grown on SrTiOa (100) viewed along the (110) azimuth, (b) SrVOa film grown on SrTiOa (100) viewed along the 
(110) azimuth, (c) Nb02 film grown on SrTiOs (111) viewed along the (1 — 10) azimuth. For LaTiOs and SrVOs the RHEED 
pattern shows four-fold symmetry with a weak 2x reconstruction for LaTiOs and no reconstruction for SrVOa. For NbOa, the 
RHEED pattern shows six-fold symmetry due to the existence of three symmetry-related rotational domains. 
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0 Is-i- V 2p spectrum of SrVOj as-grown 



(a) 


0 Is 4 V 2p spectrum of SrVOj after sputtering 
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FIG. 11: (Color Online) XPS spectrum of O Is and V 2p region of SrVOs thin film before (a) and after (b) the Ar sputtering. 
Note that even after the Ar sputtering removes surface atoms, the dP peak is still very significant. 




















